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ABSTRACT : 77i/.v paper deals with the multipoint Cauer matrix continued-fraction expansion 

(MCFE) for model reduction of linear multi-input multi-output {MI MO) systems with various 
numbers of inputs and outputs. A salient feature of the proposed MCFE approach to model 
reduction of MI MO systems with square transfer matrices is its equivalence to the matrix 
Fade approximation approach. The Cauer second form of the ordinary MC FE for a square 
transf er f unction matrix is generalized in this paper to a multipoint and nonsquat e-matrix 
version. An interesting connection of the multipoint Cauer MCFE method to the multipoint 
matrix Fade approximation method is established . Also , algorithms for obtaining the reduced- 
degree matrix-fraction descriptions and reduced-dimensional state-space models from a 
transfer function matrix via the multipoint Cauer MCFE algorithm are presented. Practical 
advantages of using the multipoint Cauer MCFE are discussed and a numerical example is 
provided to illustrate the algorithms. 

I. Introduction 

The accurate mathematical modeling or accurate identification ol linear time- 
invariant multi-input multi-output (MIMO) systems usually leads to high-degree 
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transfer-function matrices or high-dimensional state-space models. The analysis 
and design of such high-degree transfer-function matrices in the frequency domain, 
or such high-dimensional state-space models in the time domain, is not an easy task 
because of computational difficulties and implementation considerations. There is 
thus a real incentive to reduce the complexity of linear time-invariant MIMO 
models to equivalent simple ones. 

Several methods are available in the literature ( 1 - 6 ) for approximating a linear 
time-invariant MIMO system by a reduced-order model. Among them, the method 
of matrix continued fraction expansion (2 12) has the distinct advantages of 
computational simplicity and applicability to obtain both frequency-domain as 
well as time-domain reduced-order models. Chen ( 2 ) used a matrix Cauer second 
form of continued fractions to obtain a reduced-order model by retaining the first 
several significant matrix partial quotients and discarding the others. Owing to its 
equivalence to the matrix Pade approximation about v = 0, the method of using 
the Cauer second form of MCFE provides satisfactory results in the steady-state 
response only. 

To remove this drawback, Shieh and Gaudiano ( 9 , 10 ) have proposed an MCFE 
involving expansion points at s = 0 and .v = oo, which is equivalent to the matrix 
Pade approximation about s = 0 and .v = go. More recently, Chen and Hwang (7) 
have presented a multipoint MCFE method to produce better reduced-order 
models for linear time-invariant MIMO systems. 

It is worth noting that the mentioned MCFE methods are derived for the MIMO 
systems having an equal number of inputs and outputs while being described by 
their matrix-fraction description (MFD). For broadening the application scope of 
the MCFE method to include the nonsquare MIMO systems, Shieh et al. ( 11 , 13 ) 
have proposed the use of a matrix pseudo-inverse for the Cauer second form of 
the MCFE. Their method is only applicable to an MIMO system described by its 
MED and the connection of the MCFE method to the matrix Pade approximation 
method has not been explored yet. 

The main purpose of this paper is two-fold. Firstly we wish to generalize the 
CFE algorithm of multipoint Pade fitting ( 14 ) to the matrix version for obtaining 
the reduced-degree MFDs and reduced-dimensional state-space models of a non- 
square linear MIMO system characterized by its transfer-function matrix or its 
MFD. Secondly we want to reveal the connection of the multipoint Cauer MCFE 
method to the multipoint matrix Pade approximation method. 

The paper is organized as follows. In Section II, computational algorithms 
involving a matrix pseudo-inverse are derived for determining the multipoint Cauer 
MCFE of a nonsquare linear time-invariant MIMO system from its transfer- 
function matrix and MFD, respectively. In Section III, an efficient recursive algo- 
rithm is presented for obtaining a set of MFDs w'ith different reduced degrees from 
the multipoint Cauer MCFE algorithm. Also, the connection of the multipoint 
Cauer MCFE method involving the matrix pseudo-inverse to the multipoint matrix 
Pade approximation method is established. In Section IV, a canonical block state- 
realization of an MFD based on its multipoint Cauer MCFE representation is 
presented for obtaining reduced-dimensional state-space models. An explicit 
matrix-relationship between the multipoint Cauer MCFE canonical form and the 
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block companion form is also derived. Such a matrix-relationship is useful in 
connection with the construction of an aggregation matrix which relates the state 
vectors of the system and its reduced-dimensional models. 

Throughout this paper, the following notations are used : 

A * : left pseudo-inverse of the matrix A, defined as (A 1 A) ' A 1 ; 

A': right pseudo-inverse of the matrix A, defined as A' (AA 1 ) ; 

|„ : identity matrix of order n ; 

O,, : null matrix of order n ; 

0 /lx „: /j by </ null matrix; 

set of p x q real matrices. 


II. Multipoint Cauer MCFE of a Nonsquare Transfer-function Matrix 

A linear time-invariant (/-input p-output system may be described by the p x q 
transfer-function matrix 


G(.v) = (=1.2 p ; ./= 1.2 q 

or by the matrix-fraction description (MFD) ( 15 ) ; 


C,(s) = 


Z A 2.A' 


/= n 


N 

Z A,,v' 

.t 0 J 


^ A 2 (.v)A,(.v) P^q 


(I) 


(2) 


where A , ,e /?"*" and A : ,e Let the Taylor series expansions of G(.v) about 

the in distinct real points a |, <r 2 , .... <r m , be 


where 


G(.v)= Z — <T ') 7 ' ' = ’’2 m 

r-n 


I d' 

1 ° .G(.v)U„,e 
/! d.v' 


(3a) 


(3b) 


In the following, we deal with the MCFE for the nonsquare transfer-function 
matrix G(.v). Before proceeding, let p„ p 2 , .... p,„be the multiplicities corresponding 

to the distinct real points <T„<r 2 <t„„ such that p, + p 2 + ■ ■ ■ +P„, - 2jV ; These 

numbers o, will be used as expansion points taking into account multiplicities. Let 

(.V|,.v 2 , . . . ,.v 2 \) = (<r , ffi, a • • • < a m< ■ ■ ■ 


denote the ordered point system. Also, let 

v, = pi +p 2 + ' ’ ‘ +/*i i ■ 


(4) 


2. 1 . Multipoint Cauer MCFE of an MFD 

From the definition of a matrix pseudo-inverse and the MFD (2), we can write 

G(.v) as 

G(.v) = [G # (.v)] ' = [[A,(.v)A, '(•*)]*]' = [A,(.v)Af CO]'. (5) 
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Assume that Af ,) exists and let 


H, = A 1 (.v l )Af(.y l )e/r*'\ (6) 

Then, we can expand A,(s)A?(s) as 

Ai(j)Af (s) = H, +(j — i|)A 3 (5)A2 (j) (7) 

where the matrix polynomial 

i\' - i 

AT.S-) = V A Xj s\ A 3j .e/? 1x ’ (8) 

7=0 

is to be determined. Postmultiplying both sides of (7) by A,(.v) and then using the 
identity A? (.v)A 2 (i) = I,, yields the relation : 

A i (s) = H i A 2 (.r) + fs x | ) A 3 (s) . (9) 

Substituting the matrix polynomials of A,(.v), A 2 (.v) and A,(.v) into (9), and then 
equating the coefficients of like powers, we obtain 

= ^ij + 1 — H|A 2 , >+ 1 + SiA 3 ,, + ,, 7 = N— I, N—2 , . . .,0 (10) 

where A 2M = O pxi/ and A 3 ,. v = O v . 

With the substitution of (7) into (5), G(.v) can be further written as 


G(.v) - [H ! + (.v — .v , ) A, (.v) A ? (v )] ' = [H, +(.v-.v,)[A 2 (.v)A 1 '(.v)] # ]'. (11) 


Similarly, assume that A , ‘(.v 2 ) exists, and let 


H 2 = A 2 (.y 2 )A 3 '(s 2 )e R l ’ x ‘ l 

then the term in the inner brackets of ( 1 1 ) can be expanded as 


A 2 (.?)A, 1 (.v) = H : +(.v-.v 2 )A 4 (.v)A, 1 (,v) 
where the matrix polynomial 

A f 2 

A «W= I A 4 / A^e/T*" 

/- 0 


satisfies the relation : 


( 12 ) 

(13) 

(14) 


A 2 (.v) = H,A,(.v) + (.v-.v,)A 4 (.v). 

From this we have 


(15) 


A^j A 2,2+ 1 Hi A 

where A 4 ,, = O px<l for j > N— 
now expanded as 


xj+ i +.v 2 A 4 ,, + i , j= N-2,N~ (16) 

1. With the substitution of (13) into (II), G(.s) is 


Gfr) = [H, +(.r-.r l )[H : + (.v-.v 2 )A 4 (i)A 3 '(.v)] # ]'. (17) 

By repeating the expansion procedure of (5)-(17) in the MFD, A 4 (.v)A 3 '(.y), we 
can finally expand G(j) into the following multipoint Cauer matrix continued 
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fractions : 

G(A') = [H | + (S — A' ! ) [H 2 + (s — S 2 ) [H 3 + (s — A’ 3 ) [. . . [H 2jV I 

+ (.V- .v.w (18a) 

The matrix partial quotients H,can be obtained via the recurrence relations: 


fA,(v,)A* 

/= 1,3, ...,2N-\ 

(18b) 

H = < 

lA,(.v,)A, + 

‘i ( v,). / = 2, 4, . . . , IN 

;V 

A, +2 (.v)= I 
./- o 

1 

2j s> 


1 

S — Sj 

(A,(,v)-H,A (+i (a)) 

(18c) 

A • == A- - 
+ 2,/ 1 

|— H, A /4 i j + | ~\~ 2 m j r h 



/ = JV-/ 2 -l,JV-/ : -2 0 (18d) 


where /, denotes the integer part of //2 and the required matrix inverse and matrix 
pseudo-inverse are assumed to exist. For convenience, we can construct an expan- 
sion array as shown in Table I for evaluating A,,,. In this array, the first two matrix 
rows are taken from the matrix polynomials A , (.v) and A 2 (.v), respectively. The 
subsequent matrix rows are evaluated using (I8a-c). 

2.2. Multipoint Center MCFE based on Taylor-series expansions 

In this subsection, we derive an alternative tabular algorithm for obtaining the 
matrix partial quotients H, ol the multipoint Cauer MCFE (18) ol the matrix 
transfer function G(.v) from its Taylor-series coefficients T M in (3) for 

/ = 0, 1 Hi — I and / = 1 m. 

To begin with, we define 

GJs) = G a « + G a ,(.v-iT,) + G 1 -j.j(j-ff I -) 2 + ••• 

^ I <T,y m 

k - 0 
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G - 1 M = 1/ + O t/ (*v - <*,-) + (.V - a,) 2 + • ' • 

X. 

- I / = 1,2 m (20) 

A - 0 

and 

G/.2 (*v) = T, o + T , , (y — a,) + T ; 2 (.v — <t, ) 2 + ■ • ■ 

- Z G - 2.* - <*,-)*. / = 1,2,..., m. (2 1 ) 

A = 0 

Then, from (3a), (20) and (21), we can represent G(.v) by one of the following 
infinite-term MFDs: 


G(.v) = G,i(.v)G, i '(.v) 

= [G/.i(.v)G,f 2 (.v)]', / = 1, 2 m. (22) 

According to the form of the multipoint Cauer MCFE (18) and the expansion 
ordered point system (4), we know that matrix partial quotients H, + , , H„ + 2 , . . . , 
«,♦* are evaluated from the infinite-term MFD G fi2 (i)G,; ,'(.?). Since the multi- 
point Cauer MCFE involves m distinct expansion points cr„ i = 1,2, m, it is 

therefore required to manipulate the m infinite-term MFDs in (22), which represent 
the same transfer-function matrix G(.s), to obtain the partial quotients H, of the 
multipoint Cauer MCFE (18). 

Now, starting from (22), we expand the infinite-term MFD G, 1 (.v)G, # 2 (.v) as 

G,,(s)G* 2 (.v) = H , +(.9-.v,)G,,(.v)G5(.v), / = 1,2 m (23) 

where 

Hi = G u CV| )Gf 2 (.S'i ) for .v , =rr, 

~ Gi t i >0 Gf 2,0. (24) 

In (23), each of the infinite-term series C /f3 (.y) satisfies 


G/. iCv) — H|G ( 2 (.v) + (.v — ,V|)G, ,(.v), i~ 1,2 m (25) 

which is obtained by postmultiplying both sides of (23) with G, 2 (,s) and using the 
identity G,f 2 (.v)G, 2 (*v) = \ if . Substituting the series representations of G, ,(.y), G, ,(.y) 
and G, 3 (.y) into (25), we have 

X. / 

Z G,,.a.v-ff,)* = H, Z G, 2 ,*(i-(r,)‘ 

* = <> A O 

/ 

+ [(•'- cr, ) + (<r,-.v i )] Z G, , A (.v — ff,) A . (26) 

A - 0 

Equating the coefficients of like terms ( v — cr,)*, we obtain 
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G|.} 


[G,. i ,-H i G, 2 ,-G u , 

< a, — s , 

-G,, u + 1 — H | G, 2 .k * i - 


/ = 2, 3 m:k = 0, 1 , 2, . . . 

/= 1 ;A = 0, 1,2,... 

(27) 


where G, 3 i = 0, r 

In general, the MFDs, G,,(.v)G* . ,(.v) for /odd and G, ; (.v)G„! ,(.v) for./ even, can 
respectively be expanded as 


Gjj(s)Gfj i i(.v) = H ; + (.v — a 7 )G, /f ^ [(.v)^ 

/ = /,/+ 1 ,.. 

,,w (28a) 

and 



G,,(.v)G,,! ,(.v) = H / +(.r-.v / )G /Jl 2 (.r)G,. ; | , CO. 

/ = /,/+ 1 . 

,,/h (28b) 

where .v ; = a { and 



(Gij(ot)Gi j » i(rr/) = G/ (,G/ /+ i,o 

[G // ((t / )G / i/ | 1 (<r,) = G /; 0 G ul i.o 

for j odd 
for /even. 

(29) 


In (28), each infinite-term series G, ; . 2 (.v) satisfies 

G, ,(.v) = H,G,, ,(.v) + (.v-.v,)G, (+2 (.s), i = /,/+! m (30a) 


t G ilk (s-°i) k = H, I G,„ u(-»-ff/)* 

AO * 0 

+ [(.v — ff/) + (ffi X G,. ; , :.a(.v — <r,) A . (30b) 

A - (I 

This last equation allows one to evaluate the coefficients G, ;t 2 .*for/ = /, /+!,...,/» 
and k = 0, 1,2 as follows: 


[G,,.*-H,G M . u -■ (>, ,* 2j< : 1- 


G,,/-. 


.* = 1 


i = /+ l,/+2, m 


(31) 


Gw* (. i — H,G ( / f i , a ( i , i — /. 


Following the above derivation, we can construct m expansion arrays as shown 
in Table II for calculating the matrix partial quotients H,,H,, . . . , H 2iV of the 
multipoint Cauer MCFE of G(.v) from the matrix coefficients of the Taylor series 
expansions in (3). 

In the expansion array in Table II, the first two matrix rows are specified by 


fl,, for k = 0 

G ' u = {o„ for A =1,2 — I 

G, 2 .*=T,a for A = 0, 1 1 ■ 


(32a) 

(32b) 
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Tabu: II 



Expansion array for evaluating Gj 

J.k 


G„,.o 

C, M 

G,.,.2 

G, i,. , 

c„. 

C(,2.tt 

G,,.i 

G, > 

Go, , 

G l2fi 

G„.„ 

G, i.i 

G,.,, 

Go,! : 

Go,; 

| / t l 0 

\ , / 

‘v, \ 2 \ p 

! ,i,o 

G,., )t t,i 

g, v ,. 3 

G,., (t ,,^.. 2 

G ( , v 

Gj.,., | 2. , 

G;.,. , 2.2 

G,, s * 2.p, 2 

G,' 4 ; 

G,,r.,, 

G,, . 

G/.v, 4- Ijj, - 2 


I / G/. v< + ft, 1.0 

V I 1 A 1 ^ P 

| / 

G,., jf l. t 

G,>. r J ( j , i 

C 4 . l(ft ,. 2 








The subsequent matrix rows of the expansion array are obtained by the following 
recurrence formulas, 

- [G/.y.*— G,. /+ 2,* i]. 


G 


Lj + 2.A 


< A = 0, 1 /i f — 1 and /' = 1 , 2, . . . , v, 

1 1 H ( G, . |X' f i, 

A- = 0, 1 v, + n, -j— 1 and j = v, + 1 , v, + 2 v, + //,— ! 


(33) 

which are derived from the relation in (30). The matrix partial quotients H, for 
j- v,+ l, v, + 2, associated with the expansion array in Table II are 
evaluated as 


fG,.,.oG* + i.o for / odd 
1 lG W0 Gyi 1,0 for j even. 

It is noted that in the foregoing derivations, we have assumed that the required 
matrix inverse and matrix pseudo-inverse exist. If G,J, or G,% for i= 1,2,..., m 
do not exist, we should readjust the expansion points o , such that the matrix inverse 
exists. When all coefficients of any matrix row in the expansion array are null, the 
expansion procedure terminates. In this case, we can use the calculated matrix 
partial quotients and the inversion procedure in the next section to derive the M FD 
for G(.v). 


///. Reduced-degree MFDs Derived from the Multipoint Cauer MCFE 

Suppose the transfer-function matrix G(.v) of a high-degree (/-input /r-output 
system is expanded into a multipoint Cauer MCFE in ( 1 8). An A/th-degree reduced 
MFD, G w (.v), for the system can be obtained by truncating the matrix partial 
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quotients H, Wt „ H 2 M+ 2 ,..,H 2V in (18). Hence G A ,(.r) is given by 

G, ,/(.?) = [H ! + (.v — .v j ) [H 2 + (.v — -v 2 )[H , + (s — .Vj)[. . . [H,„ i 

+ (s-s 2M ,)H? A , ]'...]*] ']*]'• (35) 

To convert the truncated Cauer MCFE into its MFD, we define 


P*(.v)Q* '(.v) = < 


[H, + (.y-J l )[H 2 +( 5 -.v 2 )[H, + (.v-,v,)[...[H* , 
+ (s-s k i)Hf]'. . .] # ] '] # ] ' forAeven 
[H , + (s - s , ) [H , + (.v - s 2 ) [H 3 + (.r - ,v,)[. • -[Hi , 


+ (s-s k ,)Hr] # ...]'] # ]'] # for A odd 


(36) 


where 

p*(j) = P*.o + P*.i-«+ •■•+***.!(* n/ 2 | >ylt * U/2, < W,eR pyq (37a) 

Q*(.v) = Q*.o + Q*.i- y + ‘ ‘ ‘ +Qm* Qu e R " ‘‘ (37b) 

and [r] denotes the integer part of r. Then it can be shown that the matrix poly- 
nomials P t (.v) and Qt(.v) satisfy the following recurrence relations: 

P A (.v) = P* | (.v)H* + (.v — s k ,)P* 2 (.v), A =1,2,... (38a) 

Q. (.V) = Q k , (.v)H* + (s-s k , )Q k 2 (s) , A = 1,2,... (38b) 

where P„(.y) = O,*,, P,(.v) = O,,. Q„(.y) = I„, and Q,(.y) = H, Using the relations 
in (38). we can construct an inversion array as shown in Table III to compute the 
matrix coefficients P kJ and Q u from the matrix partial quotients H ; . 

In the inversion array, the first two matrix rows are given by 

Po.o = <w P,.o = I,.< Qo.o = Qt.o = H, (39) 


TABU: III 
Inversion array 


Numerator 


Denominator 


P = 

P d.n = 

P,.o 


O p 

o„ 


P,.o 

P 4 1 1 


P., 

P4. 


P A 


P* 


IH 


Q 1.0 F).y + / . 

Qo.o = 1 „ 

Q1.0 

Q:.o Q2, 

Qjii Qu 

Q4.0 Qu Q. 

Quo Q1.1 ••• P*.i* a 
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and the remaining matrix entries are evaluated using the following recurrence 
relations: 


P M = P * ■ ijH*— .** ,P* 2,/ + P* 2 .j 1 1 7 = 0,1 [(A — l)/2] (40a) 

Q*/ = Q* iQ* 2j + Q< I., i > 7 = 0,1 [A/2] (40b) 

where 


and 



'0/.X, 

for k even 

P, 

for k odd 


Q* 


Q ( , for A- even 
O i; , p for A odd 


for A = 1 , 2, . . . , M. 

Once the inversion array is developed, the A/th-degree reduced MFD, which 
corresponds to the truncated MCFE in (35), can be readily written as 

C„(.v) = F* 2 ,/ (.v ) Q 2 / (v) 


= [Pl«.0 + P 2 ,v/, 1 .V+ • • ■ + P 2 A/ . M ,.V A ' ' ] [Q 2A/.0 

+ Q2M,l '’+ ‘ ' ' +Q 2 M,W^ W ] (41) 


It is important here to examine the link between the A/th-degree reduced MFD. 
G m (a) and the original transfer-function matrix G(.v). As a specific example, we 
consider G,(.v) which has the following explicit form : 


6|(.y) = H,[H|H 2 + (.y-.v l )IJ (42) 

With H, being represented by (34), we have the following explicit expressions for 
G t (,v | ) , G,(,v 2 ) and (d/ds)G ,(.v , ) : 


Gj(.v,) = H : [H,H ; ] 1 


I.2.0G I 3 o[G U,oGf 2.0 G K2.oGi.io] 

1 for ,y 2 = s 

5 . 2, 0 G 2 ],() [G M.oGf 2 ,oG 2 ( 2 ,()G : jo] 

1 for s 2 / .v 

r 2 ) = H 2 [H 1 H 2 + (.v 2 -.v l )IJ 1 


| G 2 , 2,0 (*v 2 — *v i )G 3 

*.«>] ' 

= G 2, 2 ,()G 2, ],o — G 2 , 2.0 


H 2 [(H,H 2 ) 2 ] 1 


G|. 2 .()G L3 „ = G j 2 ,(|G | 2 .(|G | 2, j 

for s 2 = s ] 

~ G2.2.0G2. .],(! = [(Gf i 2 ,«C 2 .2.o)"] 1 

for .v 2 ^ .y, 


(43a) 


(43b) 


(43c) 
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Comparing (43) and (21), we have 



G,(.v,) = G(.v, ), G,(.r 2 ) = G(.v 2 ) for 

,v j 4- s 2 and p = q 

(44a) 

G,(.v,) = G(.v l ), ^G,(.v 2 ) = d d v G (.V:) 

for s , = s 2 and p = q 

(44b) 

G> , (,v , ) # G(.v l ), G,(.v 2 ) = G(.v 2 ) for 

P # <■{■ 

(44c) 


By a similar direct deduction, we can obtain the link between the A/th-degree 
( M < N) reduced MFD, G„(.v), in (41) and the original transfer-function matrix 
G(.v) as follows : 

(i) square system (p — q) 

d ' G m (.v)|, „ = * G(.v)|,-.„, for ./ = 0, 1,-- ..fi,- I and /= 1,2 / 

d.v' ' av 

(45a) 

where / is the number of distinct points that appear in the ordered point system. 


and 

(.V,,A’ 2 S 2 .u) = (<T, (7/ ,, (T f <T, ) 

ft , )/, ft : -- f* : to ' -to •. to ' 2 M v, 

(45b) 


which is obtained by truncating the last 2N — 2M points of the ordered point 
system in (4) ; 


(ii) nonsquare system (/) ^ q) 


d w G.v/(.v)l, = dy/ GCv)l. v =ff, for fij > 0 , j = 0, 1 fir 


1 and 1 ^ / < / 


(46a) 


vhere fi r the multiplicity of the point <x„ appears in the ordered point system. 


(.V 2 , .V 4 , S 2M ) = ((T\ ^2 (T h ...,(Ti 

a to 


(46b) 


which is obtained from (45) by deleting the odd-indexed elements .v 2 , , for 
/= 1,2 M. 

Notice from (46) that for a nonsquare system, the A/th-degree reduced MFD 
obtained by the multipoint Cauer MCFF. has the same Pade properties only around 
the even expansion points. However, since the multipoint Cauer MCFE involves 
the matrix pseudo-inverse in the calculation of odd matrix partial quotients, the 
obtained reduced-degree MFDs approximate the Pade properties around the odd 
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expansion points, i.e. 

d' * d' 

^G M (i)|„ ti » d?/ G(,v)L. iai for y = 0 ,l, and 1 </</ ( 47 a ) 

where the expansion point a, appears fi,- — /Z,-( > 0 ) times in the ordered point system, 
and 

(.V, ,.V 2 , . . . ,S 2M l) = (^l, • * • , G\* 0*2, . . . ,<J 2 , . . . , <T h . . . , G { ) (47b) 

L V j ' ¥ * X v V 

th fi: ft/ ji t 

which is obtained from (43) by dropping the even-indexed elements ,sv for 

7 = 1,2 M. 

IV . State-space Models Obtained from the Multipoint Cauer MCFE 

4.1. State-space formulation of the multipoint Cauer MCFE 

Let Y(.s) and U(.v) be, respectively, the Laplace transforms of the output and 
input vectors of the system. Then Y(.v) is related to U(.v) by 

Y(.v) = G(.v)U(.v). (48) 

Premultiplying both sides of (48) by G # (.y) and using (18), we have 
U(.v) = G * (.v)Y(.v) 

= [H,+(.v-.v l )[H : + (.v- .v 2 )G i (v)]' > ']Y(.v) 

T 1 “| # 

- H,Y(.y)+ ^ ^ (H, +r l G,(v)) + G,(v) I Y(.v) ( 49 ) 

where 

C s ) — H / + ^ ^ y 1 r ! ^* ,+ 2 ( >v ) ) “h ^i + : ( s ) 

and r, = >v,-.v l+ , for / = 1 , 3, 5, ... , 2 N- 1 . Let 

T 1 

^i(‘ v ) = v _ y (H 2 + r | G ^ (.v) ) + G 3 (.v) | Y(\). (50) 

Then, we have 

Y(*v) = 1 (H 2 +r 1 G,(.v)) + G,fv) X t Cv) (51a) 

- — * s i 

and 

X i (.v) = U(.v) — H i Y (.v). (51b) 

The above two relations allow one to construct the block diagram as shown in Fig. 

1 for representing the relationship between Y(x) and U(.v). 
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Fig. 1. Block diagram representation of G ,(.?). 


Using a similar procedure to expand G 3 (.v), Gs(.v), . . . , a realization block diag- 
ram (Fig. 2) for G (s) can be achieved. It is readily apparent from Fig. 2 that a 
canonical state-space representation of G(.v) can be written as follows: 

z = Hz + Gu (52a) 

y = Fz (52b) 


with 


z = 

[*T 

Z 2 

... Z T vl T 

, z, 

H = 

— H (I H C + R ® I„ 




H, 



Ho 

= 

H, 

H 3 

O 



H, 

H, 

H2A 1 



r^! 






r? 

O 

R 

= 

S\ 

•v 3 r, 




S\ 

S 5 

r IN 



f l 

S 3 5 s 

. S 2 , V 

G 

= 11, 

I, ••• I, 

]’ 



A' block elements 


F 

= [H 2 

h 4 

H 2 ,v] 




*J T 


-1 

h 2 

h 4 

Mv 

, H e = 

0 

h 4 





H«. 


3 

3 f 2S I 


where the superscript T designates the transpose of a matrix or vector, and ® 
denotes the Kronecker product ( 16 ). 


4.2. Link between MCFE and block companion forms 

Let A, v = 1„, then a controllable block companion form for the realization of 
the MFD (2) can be written as follows ( 17 ) : 


Vol. 331 B, No. 2, pp. 189-216, 1994 
Printed in Great Britain. All rights reserved 


201 






Tong- Yi Guo et al. 



Fig. 2. Realization block diagram of multipoint Cauer matrix CFE. 


where 


x = [x 


A = 



x = 

Ax + Bu 



y = 

Cx 


T x] 

... xJ] T , 

x, = k. 

X it 2 

O, 

h 



o„ 

Ov 

L 

O 

O, 

o, 


K 

— A 1.0 

— A | | 

"A 12 

~~ A 


... ij T 




. . x iA 


C — [A 2,0 A 2 


(53a) 

(53b) 


Assume that the system is completely controllable with controllability index equal 
to Nq. Then there exists a nonsingular block transformation matrix P which links 
the block state vectors z and x by 


z = Px (54) 

and satisfies the following relations: 

PA = HP (55a) 

PB = G (55b) 

G = FP. (55c) 

The desired block similarity transformation matrix P can be derived from the 
relation (18a) and (55c). Let the block similarity transformation matrix P be of 
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the form 


p.,1 

P.2 

P \.N 

p 2 . 

P 22 . 

P 2.N 

Pv 1 

P.V.2 

P.\’,;V 


(56) 


where the block elements P,., are q x q matrices. Define the matrix polynomials 

P,(s) = X P ..r v> ' ’ /' = 1 , 2, . . . , (V. (57) 

j= i 

Then the relation (55c) is equivalent to 

A 2 (x) = H 2 Pi(.9) + H 4 P 2 (.v) + ■ • +H2,vP,v(- ? )- (^8) 

Rewrite the recursion in (18b) as 

A ,(.v) = (.9 — .s,) A,+ 2 (x) + H,A, + , (. 9 ), / = 1,2,.. . ,2 N. (59) 

Then, by applying the above relations recursively for i = 2,4 2N and noting 

that A 2 , v+2 (s) = O pxll , we have 

A 2 (.v) = H 2 AjCv) + (s-S2)H 4 A 5 (.y) + (j-s 2 )(5-54)H 6 A 7 (s) + ••• 

+ {S — S 2 ) ' ' ' {•'>— s 2N 2)H 2 jvA 2(v+ ,(.?). (60) 

Comparing (58) and (60), we obtain 

p ( (, 5 ) = 1] (j-j 2y )A 2<+ ,(j), i= \,2,...,N. (61) 

j= I 

Using (57) and equating the coefficients of the like powers of s' in (61), and using 
the relations in (57), we have the explicit form for P, v as 

P u = X n,k A 2 i+ \j-k+i, / = 1, 2, . . . , n and y = 1, 2, . . . , « (62a) 

k= \ 

where A,,, are the coefficients of the matrix polynomials defined in (18b,c), and r, M 
are given by 

! 0 for i < k or k = 0 

1 for i = k (62b) 

r,. M ,-.v, 2 r t i,k for i > k. 


4.3. Reduced-dimensional state-space models 

Truncation of the last (2N-2M) matrix partial quotients of the multipoint 
Cauer MCFE in (18) after the first 2 M matrix partial quotients is equivalent to 
discarding the N-M inner block loops of Fig. 2. Hence, the multipoint Cauer 
MCFE canonical state-space model corresponding to the Mth-degree G, w (.v) can 
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be readily obtained as 

z = Hz -f Gu 
y = Fz 

where 


(63a) 

(63b) 


Z=[zj zl ... H]\ Z, = 

H = — H 0 H e + R (x) \ q 


H 

1 


H 

h 3 

O 

H 

h 3 

H 2M- 1 _ 

~*i 



**1 

ry 

O 


Sy r 5 


S\ 

*3 ^5 

r 2M 3 

Si 

*3 ^5 

S 2M- 3 


I, ... 

y T 


M block elements 


F=[H 2 h 4 H 2v/ ] 


1 ^i,2 ••• 




h 2 

h 4 

h 2M 

a = 

o 

h 4 

h 2m 




h 2 „ 


r 2 M-i 


It is obvious that H is the upper-left M by M block submatrix of H, and G and F 
are the upper block subvectors of G and F, respectively. Hence, the truncation 
procedure in the time domain can be represented by the following mathematical 
relations : 

H = JHJ 1 (64a) 

G = JG (64b) 

F = FJ T (64c) 

where 


“ Whlq : M)q\- 

Based on the multi-feedback multi- feed forward control theory ( 2 , 10 ) and the 
realization block diagram of Fig. 2, we know that the inner block states have fewer 
contributions to the system output. Hence, it is reasonable to consider that for an 
appropriate A/, the following approximation holds 


z % Jz. (65) 

With this approximation and the relation (64), we conclude that the state vector z 
of the reduced-dimensional state-space model is approximately related to that of 
the original system in (53) by 
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z = Qx (66a) 

where 

Q = JP. (66b) 

The Q is called an approximate aggregation matrix. The availability of this 
approximate aggregation matrix enables the reduced-dimensional state-space 
models obtained by the multipoint Cauer MCFE method to be useful in designing 
reduced-order controllers. 


V. An Illustrative Example 

To illustrate the proposed algorithms, we consider a two-input-three-output 
system characterized by its third-degree MFD 


- 

- -7.6051 

-18.4986' 


"-23.5 164 

-49.8596“ 


-35.0581 

-85.4861 

+ 

-4.5654 

21.8066 


_ — 25.4654 

— 62.0624_ 


-6.4341 

8.5709_ 


r -6.6984 


+ 


1 .8038 - 

0.2685 


" -64.1796 
-167.1634 


3.8209“ 

0.8574 



- - 

0.0243. 



L L 


- 148.9154“ 
-326.1264 


- 1.2972 -2.5216“! 

-85.9888 -209.61 loj 

“- 16.0738 29.8019“| 

_ 39.1889 40.4941 J 


2 


+ 


1 

0 


0 

1 



I 


The transfer-function matrix is given by : 

r<7 u Cv) 0i, 2 (O 


G ( 5 ) = 


02.1 CO 02,2(0 

. 03 . 1 (0 0 3 , 2(0 J 


where 

9u (s) = (3.440200 + 29. 863961.V+ 34.330077r -90.583719.V 4 

— 145.02521 l.v 4 — 6.6984.v 5 )/c/(.v) 

g 2 , (_*,-) = (- 2.283757 -24.693759.v- 87.4 12685.V 2 - 96.945897s 1 

+ 34.877295.V 4 + 1.8038s 5 )/</(.v) 

g }J (s) = (1.156658+ 1 6.036971. v + 1 3. 523988.V 2 — 33.6272 1 4.v 3 

+ 5.390856.V 4 + 0.2685 j 5 )/r/(.v) 
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9\A*) = (4.8 19364 + 60.095 159s + 200.1 57 141s 2 + 241. 04929 Is 5 

+ 88.349065s 4 + 3.8209s 5 )/<f(s) 

9iA*) = (22.490064 + 225.973085^ + 345. I479I3.v 2 +23.6977555' ; ' 

— 18.168391s 4 — 0.8574s s )/rf(s) 

(hAs) = (1 6.293793 + 163.607382.? + 248. 927928.V 2 + 30.343293a- 5 


+ 0. 1 78497s + 0.0243 s 5 )/t/(s) 

d{s) = 55.07803 l' + 649.225525s+ 1817.984559s 2 + 1578.243399s 3 

+ 1 26.7036 1 4s 4 + 24.4203s 5 + s 6 . 

The controllable block companion form corresponding to this MFD is 


[i::] 


"0 0 " 
_0 0 


n 0-1 

Lo ij 

“0 0 “ 
_0 0 _ 


r° °i 

_o oj 


1.2972 2.5216 

85.9888 209.6110 


;][ 
:] 
:] 


64.1796 148.9154 

167.1634 326.1264 


:] 


16.0738 

39.1889 


-29.8019 

-40.4941 


'] 


t:;] 

C:;J 


+ 


[::] 

[::] 

[::] 


[:;] 


yr 



- -7.6051 

- 18.4986“ 

y 2 

— 


-35.0581 

-85.4861 

~y 3 - 


_ 

.-25.4654 

-62.0624, 


-23.5164 
-4.5654 
L -6.4341 


-49.8596 

21.8066 

8.5709 


“-6.6984 

3.8209“ 

- 

1.8038 

-0.8574 


. 0.2685 

0.0243, 





The Taylor series expansions of G(s) around expansion points at s = 0, 1,2 are 
evaluated as 
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~ 0.062460 0.08750 11 r- 0.1 94033 0.059688" 

G(.s-) = -0.041464 0.408331 + 0.040411 -0.710370.? 

0.021000 0.295831 J L 0.043629 -0.516607. 

' 0.848784 0.042327“ 

+ -0.694785 1.161968 ,? 2 + • ■ • 

_ -0.961 898 0.844367. 

0.04 1074 0.1406861 r- 0.09 1768 0.046501“ 

= -0.041070 0.140685 + 0.012861 -0.105319 (.v — I ) 

0.000647 0.108021 J L- 0.009672 -0.068097. 

0.005283 -0.005326“ 

+ 0.008434 0.058837 (v— l) 2 + * - 

.0.008503 0.042513. 

--0. 126800 0.1819711 [“-0.078945 0.036070“ 

= -0.023260 0.071566 + 0.020381 -0.046636 (,?-2) 

.-0.003600 0.065727 J L- 0.000845 -0.026606. 

"0.007431 -0.005244“ 

+ 0.000944 0.014180 (i-2) 2 + ■ ■ ■. 

.0.002074 0.009571. 

Choosing the ordered expansion point system as (s ,, s ,, .? 2 , s 2 ,s 3 , .? 3 ) = (0, 0, 1 , 1,2, 2) 
and using the recursive algorithms ( 1 8) and (33), we construct the expansion arrays 
as shown in Tables IV and V for evaluating the multipoint Cauer MCFE from the 
MFD and Taylor series expansion of G(.v), respectively. 

We obtain from (18a) and (34), the same set matrix of partial quotients as 

follows : 





0.029679 

0.095386“ 

“10.781239 

-5.586419 

4.521980“] 

-0.019033 

0.224562 


H = 

0.550334 

1 .447022 

1. 220228 J 

0.010361 

0.196506. 






-6.737767 

3.651197- 

'- 10.576133 

-5.571830 

-4.02142 1“| 

h 4 = 

1.817531 

- 1 .093386 

4.017897 

-8.535413 

- 4.024899 J 


0.235814 

— 0.184328. 






"0.009688 

0.074317“ 

“165.364859 

70.422714 

271.800239“j 

h 6 = 

0.005302 

0.01 1424 

209.425455 

49.768519 

1 29.772792 J 


_0.022325 

0.012122. 
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Table IV 

Expansion array for evaluating Aj j 


-64.1796 

— 148.9154 

16.0738 

29.8019] (" 1 o" 

-167.1634 

-326.1264 

— 39.1889 

40.4941 J 10 1 

-23.5164 

-49.85961 

r- 6.6984 

3.8209“ 

-4.5654 

21.8066 

^ 1.8038 

-0.8574 

-6.4341 

8.5709J 

_ 0.2685 

0.0243_ 

65.005884 

-16.291817] 

~1 o' 


-38.440311 

39.602353] 

0 1 


-6.728079 

3.725514- 



1.822833 

-1.081962 



0.258139 

— 0.172206_ 




"l 0" 
0 1 
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Table Va 

Expansion array for evaluating G, j k 


0 0 
0 0 


I 0 
0 I 

- 0.062460 0.08750 T 

-0.041464 0.408331 

0.021000 0.29583 1_ 

2.120382 -2.275849 

-0.004929 1.625451 


-- 0.1 94033 0.059688 

0.040411 -0.710370 

0.043629 -0.516607. 


Table Vb 

Expansion array for evaluating G 2 j.k 


1 0 
0 1 _ 

- -0.041074 0.140686' 

-0.041070 0.140685 

0.000647 0.108021. 

"1.210471 -1.219319" 

0.081244 0.587192 

0.084749 0.120865- 

-0.036275 -0.014384 

.-0.027860 0.005267. 

0.308025 -0.108965" 

0.545205 —0.167977 


[So] 

-- 0.09 1768 

0.046501 ~ 

0.012861 

-0.105319 

_ — 0.009672 

— 0.068097- 

' — 0. 1055 14 

0.4375591 

_ — 0.037548 

-0. 377289 J 

- -0.000306 

-0.051362" 

0.055560 

0.002118 

0.026659 

— 0.003759_ 
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Table Vc 

Expansion array for evaluating G 3 j k 


1 0 

0 1 



0.1 26800 

0 . 181971 " 


-- 0.078945 

0 . 036070 " 

- 0.023260 

0.071566 


0.020381 

- 0.046636 

.- 0.003600 

0 . 065727 _ 


.- 0.000845 

- 0 . 026606 . 

1.1 26698 - 

- 0 . 929643 " 


["- 0.078950 

0 . 200273 “ 

0.053916 

0.358048 


[- 0.019465 

— 0 . 138975 _ 

- - 0.082691 

0.087704 


0.003973 - 

- 0 . 022161 " 


- 0.006962 

- 0.013266 


0.015105 

0.000825 


.- 0.012934 

0 . 002500 . 


. 0.008366 - 

- 0 . 001936 . 


" 0.161345 - 

0.065929 


0.080466 

0 . 028638 " 


0.274680 - 

0.097505 


— 0.147505 

0.046821 


" 0.001500 

- 0 . 000500 " 


" — 0.001 1 19 

0 . 000338 " 

0.000121 

- 0.000048 


- 0.000046 

0.000017 

.- 0.000350 

0 . 000074 . 


_ 0.000502 

— 0 . 000133 . 

0.028680 

0 . 007509 “ 




0.023921 

- 0.007708 





Once the matrix partial quotients H, are evaluated, the canonical block state- 
space model corresponding to the multipoint Cauer MCFE can be constructed 
according to (52) as follows : 


rp"T 



-0.473153 

-0.6624781 |”8 1.728620 

-44.639017“] 

UJ 



.-0.001435 

-0.61 7223 J L 0.790269 

- 0.202302 J 




0.473 153 

-0.662478“] r 22.544374 

-12.8568931 

U. J 



.-0.001435 

-0.61 7223 J L 44.324429 

-23.946844J 

p-1 



r — 0.473 1 53 

-0.6624781 [“21.544374 

-12.856893“] 

IUJJ 


- 

L -0.001435 

-0.61 7223 J L44.324429 

— 24.946844 J 


-0.175783 

-0.7922221 " 


rhuT 


T 1 oT 


-0.040246 

- 0.072222 J 


LJ 


Lo iJ 


0.046004 

0.106164“] 




r 1 01 

U 1 

0.055941 

- 0.2245 18J 


L^2,2 J 

+ 

Lo iJ 

L 

-5.997421 

-16.2825291 


K.l 


L 1 oi 


-5.134094 

— 1 5.930034 J _ 




lo iJJ 
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) ? 2 

= 


~y 3 - 


_ 


0.029679 0.095386“ 

0.019033 0.224562 

0.010361 0.196506. 


--6.737767 3.651197“ 

1.817531 - 1.093386 

0.235814 -0.184328. 


0.009688 

0.074317" 


0.005302 

0.011424 


.0.022325 

0.012122. 

_ 



Using (62) with A,, extracted from Table IV, we obtain the block similarity 
transformation matrix for relating the block state vectors z and x in (54) as follows : 


r 

192.946970 


471.696242 

[” 65.005884 

-16.2918171 

“1 o' 



_ 139.764214 

- 

340.700047. 

L- 38.4403 11 

39.602353J 

_0 1. 




"o 

o" 


T 6.043425 

16.3886941 

'l O' 




^0 

0_ 


[5.190035 

15.70551 5 J 

.0 1. 




“o 

0‘ 



'-1 


'1 o' 




_0 

0_ 



0 

-u 

_0 1. 

- 


Based on the model reduction method proposed in this paper, the second-degree 
reduced MFD model, G : (s), derived from the matrix partial quotients, H„ 
i = 1 , 2, 3, 4 is formed as 


J- 

-2.396027 1.321639“ 


-- 6.708088 3.746583“ 



- 10.902575 5.937033 

+ 

1.798498 -0.868824 

s 


-7.941521 4.336653. 


0.246175 0.012178. 

- 




-0.837195 

-26.785352 


0.6924161 -22.071221 

14.610068_r[_ — 44.322994 


13.5193711 p 0l 
24.564066J + Lo 1J _ 


£i,i(.v) .7 1 . 2 Cl 9 i.j(*) 

9iAs) g 2 ,i(s) 9i.^s). 


where 

g x = (0. 394446 + 2. 070750s — l.l 141 70s 2 — 6. 708088s 3 )/d(s) 

g xl (s) = (0.552578 + 4.730760s + 9.3 19 107s 2 + 3. 746583s , )/r/(s) 

g 2 ,(. V ) = (-0.261851 — 1.660099s— 5.233055s 2 + 1.798498s 3 )/d(s) 
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92.2(s) = (2.578663 + 1 5.8404575 + 0.798490s 2 — 0.868824s 3 )/*/(.r) 
9\\(s) = (0.132620+ 1.0602325- 1.354695i 2 + 0.246175i- , )/rf(.v) 
9x2^) = (1.868212+ 1 1.4684845 + 0.739738.v 2 +0.0121785 3 )/</(5) 
d(s) = 6.315131 +49.784100.v+70.8329265 2 + 2.4928465 3 + .s- 4 . 


Furthermore, the multipoint Cauer MCFE state-space model corresponding to 
G 2 (5) can be readily obtained as 


m-11 



— 0.473153 

-0.6624781 pi. 728620 

— 44.639017"]" 

IA 2 J 



.-0.001435 

-0.61 7223 J L 0.790269 

- 0.202302 J 

PH 



” — 0.473153 

-0.6624781 22.544374 

- 12.856893] 

. L^2,2 J _ 


- 

.-0.001435 

-0.61 7223 J |_44.324429 

- 23.946844 J_ 



~yr 



- 0.029679 0.095386- 


“-6.737767 

y 2 



-0.019033 0.224562 


1.817531 

->’ 3 - 



_ 0.010361 0.196506. 


_ 0.235814 


3.651197 

-1.093386 

-0.184328. 



The Taylor series expansions of G,(s) around .v = 0, 1 , 2 are, respectively, evaluated 
as 


' 0062460 

0.087501 - 


-- 0.1 64492 

0.05932L 

-0.041464 

0.408331 

+ 

0.063997 

-0.710663 

. 0021000 

0.295831. 


. 0.002335 

— 0.516094_ 


+ 


0.419732 

0.868085 

0.468475 


0.026594' 

1.148819 

0.867507_ 


+ ... 


0.04 1074 

0.140686" 


-- 0.09 1599 

0.046459- 

-0.041070 

0.140685 

+ 

0.012294 

-0.105179 

_ 0.000647 

0.108021. 


.-0.007988 

-0.068514. 


+ 


-0.006230 

0.008666 

.0.007286 


-0.005601* 

0.058842 

0.042725. 


( 5 - 1 ) 2 + •• • 
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-- 0.1 26037 

0.181728" 


"- 0.077841 

0.035700" 

-0.023818 

0.071745 

+ 

0.019754 

-0.046413 

0.0028 14 

0.065532. 


.-0.000689 

-0.026610. 


+ 


"0.007584 
0.000869 
_0.00 1655 


— 0.005306“ 
0.014214 
0.009714. 


(s — 2) 2 + • • •• 


Comparing the above series with those of G(.v), we can verify that the partial matrix 
Pade properties of G (s) given in (46a) and (47a) match those of G(x). 

For comparison of the time responses between reduced-degree models and G(.?), 
the second-degree reduced MFD obtained by the single-point Cauer second form 
of the MCFE ( 12 ) is given below : 


-2.446550 1.375684" 


"-6.893735 3.936407" 


-11.180131 6.225057 

+ 

1.825865 -0.899738 

.9 

-8.136202 4.539832. 


0.213142 0.045311. 

- 

-0.711690 0.584829 

-27.452343 15.304513 

(“-22.590895 14.079862 

+ L — 45.510647 25.793547. 


1 o" 


i 

9*. i (s) 

gU*) 

0*3 CO 

-0 L 

s~ 


j* iC0 

9*. 2 CO 

g*A s )- 


where 

g * ,(.*) = (0.322474 + 2.061388.9— 1.11 1976.V 2 - 6.893375.v , )/r?*(.s) 
ft 2 (s) = (0.45 1753 + 4.599323.9 + 9.511 552.v 2 + 3.936407.V 3 )/ d*(s) 
ff* iC s ) = ( — 0.21 4073 — 1 .8248239 — 5 .032249.9 2 + 1 ,825865s 3 )/ d*(s) 
g* 2 (s) = (2.108154+ 16.3576199+0.8430I3.9 2 — 0.8997389 , )/<?*(9) 

^*,( 9 ) = (0 . 1 08422 + 1 .255 1 399 — 0. 576348.V 2 + 0.2 1 3 1 429 1 )/ d (s) 
g* 2 ( s ) = (1.527333+ 1 1.8408289 + 0.5151899 2 + 0.04531l9 3 )/<i*(9) 
d*(s) = 5 . 162856 + 49. 0414889 + 72. 677l369 2 + 3.20265l9 J + .v 4 . 

The values of the impulse response energy (18) for transfer function elements and 
outputs of the original system and those of the reduced-degree MFD are listed in 
Table VI. The values of the integral squared-error of the impulse and the unit-step 
time responses for each output of the reduced models are shown in Table VII. The 
above values are usually taken as criteria for choosing proper reduced-degree 
models. The unit-step time responses for each output of the original system G(.?) 
and the reduced-degree models G 2 (.9) are plotted in Fig. 3. It is observed that the 
time-responses of the reduced model are close to those of the original system. 
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Table VI 

The values of impulse response energy of the original and reduced- 
degree models 



G*(s) 

G 2 w 

G(j) 

9 i.iCO 

9.493938 

12.707962 

14.771495 


0.748100 

1.037583 

1.283513 

9i.i(s) 

0.008769 

0.022492 

0.042638 

9 2(2) 

3.219912 

4.134660 

4.683292 

92.2 (s) 

0.298879 

0.382659 

0.448826 

92.2 (S) 

0.027614 

0.032275 

0.037869 

}'i(s) 

2.039829 

2.888713 

3.421429 

y 2 (s) 

0.191937 

0.267885 

0.331350 

?}(*) 

0.030720 

0.032574 

0.038061 


Table VII 

The integral squared-errors of the impulse and unit-step responses between the original and 

reduced-degree models 


Impulse response Unit-step response 



y i(s) 




T'jOv) 

yM 

G!(.v) 

0.467087 

0.045071 

0.003970 

0.006885 

0.000700 

0.000136 

G 2 (s) 

0.135860 

0.014271 

0.001245 

0.002030 

0.000257 

0.000098 



Fig. 3. Unit-step responses of the original and reduced-degree models. 
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VI. Conclusions 

A matrix continued-fraction methodology has been established for modeling and 
model reduction of multivariable systems having an unequal number of inputs and 
outputs. The developed methodology is based on the multipoint Cauer matrix 
continued-fraction expansion and the matrix pseudo-inverse, and is applicable to 
obtain both the frequency-domain reduced -degree MFDs and time-domain re- 
duced-dimensional state-space models. The connection of the presented multipoint 
Cauer MCFE method involving matrix pseudo-inverse to the multipoint matrix 
Pade approximation method is also established. Although the algorithms in the 
paper are derived for a right MFD, the extension to a left MFD is straightforward. 

The use of a multipoint continued fraction expansion gives good reduced-order 
models that provide satisfactory approximation of the original system not only in 
the transient response but also in the steady-state response. Although the stability 
of the reduced models by the multipoint Cauer MCFE is not completely guaranteed, 
the flexibility of choosing expansion points provides the highest probability for 
finding stable reduced-order models (19). 
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